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Abstract 

Superconductors are considered in view of applications to atom chip de- 

\ vices. The main features of magnetic traps based on superconducting wires 

' in the Meissner and mixed states are discussed. The former state may mainly 

be interesting for improved atom optics, while in the latter, cold atoms may 

provide a probe of superconductor phenomena. The properties of a mag- 
00 ' 

q . netic side guide based on a single superconducting strip wire placed in an 

external magnetic field are calculated analytically and numerically. In the 
mixed state of type II superconductors, inhomogeneous trapped magnetic 
flux, relaxation processes and noise caused by vortex motion are posing spe- 
^ ' cine challenges for atom trapping. 

PACS: 37.10.Gh Atom traps and guides 85.25.Am Superconducting device 
characterization, design, and modeling 



1 Introduction 

One of the high priority goals in atom chip research is to increase lifetime and 
coherence time for ultracold atoms trapped in magnetic potentials close to the 
surface. This is important for both scientific aims and technological applications. 
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Progress towards this goal demands the control and reduction of magnetic noise 
produced by the metallic components of the atom chip. Randomly fluctuating 
magnetic fields are generated by thermal current noise in the conducting chip ele- 
ments and reduce the number of trapped atoms (losses), increase their temperature 
(heating) and lead to a phase uncertainty in the atom's state (decoherence) - see, 
for example, 01 SO and references therein. Theoretical analysis of the magnetic 
noise generated by a normal metal JH[5l[6ll3 predicts a fast reduction of the life- 
time with the decrease of the distance z t between the trapped atom and the metal 
surface (trap height); this is in excellent agreement with lifetime measurements 
jH |H [10l|. At a trap height less than 10 20 /iin, thermal magnetic noise exceeds 
all other harmful influences on the atom cloud (technical noise due to the current 
supply instability, residual gas collisions, stray magnetic fields) and provides the 
dominant limit for the lifetime 01 El. In the last few years, the application of su- 
perconducting materials to atom chips has been widely discussed as a perspective 
to extend the lifetime of cold atoms 0[Tj3[T2l[T3l[T4]|. A recent theoretical esti- 
mate |[i~3l of the magnetic noise caused by a superconductor in the Meissner state 
showed that the lifetime of atoms trapped above a superconducting layer would 
be, at least, six orders of magnitude longer than above a normal metal in the same 
geometry. The analysis presented in [14J predicts an atom lifetime of 5000 s at a 
trap height of 1 /im. For comparison, at the same height in a normal metal trap the 
lifetime is less than 0.1 s BU. At larger heights, the lifetime in a superconducting 
niobium chip can be much larger (up to 10 11 s at temperature T = 4.2 K). 

Two first realizations of atom chips with superconducting elements have been 
reported in Refs. lfT5l and lfT6l . In both setups, the trapped atoms were 87 Rb. In 
the experiment by Nirrengarten et al. 10311 . the current-carrying wires (in "U" and 
"Z" shape) were made of niobium and operated at about 4.2 K. The obtained 
atom spin relaxation time was estimated as 115 s. This value is comparable to the 
best one achieved for atoms trapped near normal-metal wires ifTTl . In the second 
experiment lfT6l . special efforts have been undertaken to reduce the influence of 
technical noise. Using a MgB 2 film, a "Z"-shaped wire was fabricated as a part of 
a closed superconducting loop and operated in the persistent current regime. This 
permits to disconnect the current supply and get rid of its instability, i.e. techni- 
cal noise. To our knowledge, in both experiments the trap lifetimes were limited 
by processes other than the magnetic noise generated by the superconducting ele- 
ments of the atom chip. 

Aside from magnetic noise reduction (thermal and technical), the application 
of superconductors in atom chips may be advantageous due to high current den- 
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sities without Joule heating, and practically zero electric fields across the super- 
conducting elements. Atom traps with high currents and strong confinement, as 
needed for some applications, may not be operated in the Meissner state because 
of too low critical parameters. One can then use superconducting wires in the 
so-called mixed state where magnetic flux partially penetrates the chip structures. 
This is an issue that we address in this paper, complementing previous approaches 
that focus on the Meissner state [|20l . It is well known that in addition to atom 
optics, one may use cold atoms as sensitive probes of current distribution and 
noise in the nearby surface of the atom chip |[T8l [T9ll . As the mixed state of a 
superconductor exhibits vortex phenomena and flux noise much higher than that 
of the Meissner state, it may be an interesting object to study in this context, be- 
yond the advantages of carrying larger super currents. The results reported here 
identify a fairly large parameter window where atom chips can be designed with 
superconducting elements, both in the Meissner and mixed states. 

In the following Section, we review typical properties of superconducting ma- 
terials, in particular critical parameters in view of atom chip applications. Sec- 
tion [3] analyses side guide traps formed by combining a bias field with the su- 
percurrent of a rectangular wire. We show analytically and numerically how the 
current distribution is significantly modified due to screening and flux penetration 
and identify the consequences for trapping and transporting atoms on a chip. We 
also give an overview on magnetic noise in mixed-state superconductors. Sec- 
tion |4] discusses the confinement of the magnetic trap, compared to normal metal 
wires. The analytically solvable case of a cylindrical wire and the method used 
for numerical calculations are described in the Appendices. 

2 Typical superconductors 

Let us briefly survey the properties of superconductors which are important for 
building a superconducting atom chip. At low enough fields and temperatures, 
superconductors exhibit the Meissner effect: magnetic fields are expelled from 
their interior. In this regime, the field penetrates into a superconductor only over 
a small depth A from the surface - the London penetration depth. At T = 0, A is 
of the order of 50 nm and increases with temperature as [1 — (T/Tc) 4 ]" 1 / 2 , where 
T c is the critical temperature. 

Superconducting materials are classified as type I and type II superconductors 
that differ in their behavior as the magnetic field is increased. Type I supercon- 
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ductors (pure metals as Pb, Al, Hg, In) are in the Meissner state over the whole 
sub-critical ranges of temperature T < T c , external magnetic field B < B c , and 
current I < I c . Their critical parameters are quite low, however, which is the main 
problem for applications. The highest values are observed for lead: T c = 7.2 K, 
critical magnetic field B c = 0.055 T at T = K, and surface (sheet) critical cur- 
rent density J c = B c j /x = 4.4 • 10 4 A/m in zero external magnetic field at T = 
ll2i~ll . (/x is the free space permeability.) 

Table 1: Critical parameters of selected type-II superconductors. The critical current 
density values are referring to the highest quality superconducting films and tapes (Nb 3 Sn 
and (Pb,Bi)2Si"2Ca2Cu30io). Below the first critical field B c \, the material is in the 
Meissner state (vortex free), between B& and B& in the mixed state (with vortex pen- 
etration). The temperature is 4.2 K unless otherwise quoted. Data collected from Refs. 

OH M H2 M 122 M W\ M M M- 



Superconductor 


T C (K) 


B c i (mT) 


B c2 (T) 


jc (A/m 2 ) 


Nb 


9.3 


140 


0.28 


5 • 10 10 (B = 0) 


Nb 3 Sn 


18 


40 


27 


6 - 10 10 (B = IT) 


MgB 2 


39 


30 


15 


3.5 • 10 11 (B = 0) 


YBa 2 Cu 3 7 _5 
(YBCO) 


92 


25 (B\\ab, T 0) 
90 (B\\c,T ^0) 


>100 (77 K) 


7.2 • 10 11 (B = 0) 


(Pb-Bi) 2 Sr2Ca 2 Cu30io 


108 


13 (B\\ab, T 0) 


>100 (77 K) 


ps 10 10 (B = 0) 



Much higher critical parameters are observed in type II superconductors. We 
give typical examples of this type in Table [H including niobium, its compound 
Nb 3 Sn, as well as the high-temperature superconductors (HTSC) YBa 2 Cu 3 7 „5 
and (Pb-Bi) 2 Ca 2 Sr 2 Cu 3 Oi . These materials exhibit the Meissner effect below 
the lower critical field, B < B c \, which is relatively small (50 -j- 500 G at OK). 
The largest value at T = is found for niobium (B cl (0) w 170 mT) 11221 . 
With increasing temperature, the lower critical field falls down approximately as 
.B c i(0)[l — (T/T c ) 2 ], and other parameters (B c2 , j c ) also decrease towards zero as 
T — > T c . The details of these laws depend on the type of superconductor and its 
fabrication EE]. In Tabled] we present main critical parameters at T = 4.2 K. 
For highly anisotropic HTSC with layered crystal structures, the lower critical 
field (as well as the critical current density) depends on the direction of the field 
relative to the a6-planes [3l |. In magnetic fields (external or caused by a transport 
current) higher than B c \, the magnetic flux penetrates into a type II superconduc- 
tor in the form of vortices that arrange into a more or less regular flux-line lattice 
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(Abrikosov lattice) which is pinned by inhomogeneities of the material I12TTI . Each 
vortex carries one quantum of magnetic flux $ = nU/e = 2.07 ■ 10~ 15 T m 2 (e 
is the electron charge and % the Planck constant). In the simplest (isotropic) case, 
it looks like a cylindrical tube of radius ~ A in which superconducting screening 
currents circulate around a normally conducting core of radius ~ £ (superconduc- 
tor coherence length). Both A and £ depend on temperature, with A 3> £ for most 
type II superconductors. The mixed state is observed up to the upper critical field 
B C 2 where the vortex cores merge and the superconducting state is destroyed. 

If a superconducting material is used in atom chips, it is particularly important 
that the lower critical field B cl and the critical current density j c be large. The 
second critical field B c2 for most type II superconductors is typically too large 
to be a relevant limit in the magnetic fields applied in usual atom chip setups. 
Both lower critical field and critical current density are "technology dependent" 
because they are very sensitive to crystal defects (for details see Refs. ll2~9ll32~l[33ll ). 
The j c values collected in Tabled] refer to films and tapes of the best quality. For 
example, in the atom chip experiments of Refs. llT5llT6ll . the critical current density 
was 5 ■ 10 10 A/m 2 (niobium film) and 10 11 A/m 2 (MgB 2 film) respectively. 

One specific property, which distinguishes type II superconductors from both 
normal metals and superconductors in the Meissner state, is their capability to 
freeze a magnetic flux ETTl . This effect is due to pinning forces. The vortices 
can move under the action of the transport current, when the Lorentz force j$ 
is stronger than the pinning force, which can be estimated as j c $o where j c is the 
critical current density. This flux motion results in energy dissipation and induces 
a voltage drop along the superconducting element. A voltage drop of 1 /iV/cm is 
usually taken as a criterion to define the critical current. Another mechanism of 
energy dissipation predicted by P. Anderson operates at sub-critical currents and 
is connected with thermally activated jumps of vortices out of pinning centres, 
which also generate electric fields 11341 . The motion of vortices under various 
conditions has been investigated by many authors, see for example Ref. ll32ll33l . 

3 Side guide traps with superconducting wires 
3.1 Rectangular wire in the Meissner state 

We consider magnetic traps in the "side guide" configuration to illustrate the dif- 
ferences between normal metal and superconducting chips. In this section we 
analyse a trap based on wires having the form of a thin strip which is the usual 
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Figure 1: (left) Superconducting strip carrying a current and placed in an external mag- 
netic field, z = is at the top surface of the strip, (right) Illustration of the surface current 
distribution (arrows) for the strip wire. Top: screening supercurrents in an external field 
Bo, parallel to the top wire surface (this direction is called "horizontal"). The screening 
currents are particularly large at the corners, but have opposite directions on the top and 
bottom faces; the total current is zero. Bottom: transport current without an external field. 
The current density is maximal near the corners [see also Eq.(Q])]. 

shape in present day atom chips. The case of cylindrical conductors is considered 
in Appendix |Al 

3.1.1 Horizontal bias field 

Let us start with a superconducting strip wire in the Meissner state. The strip is 
infinitely long along the y-axis, with a width 2w along x and a thickness d along 
z, perpendicular to the chip surface (Fig{T|). The strip is placed in an external 
bias field B parallel to the x-axis. If we assume that the strip is thin, d <^ w, 
the real current density distribution can be replaced by a sheet current J{x) (in 
A/m) determined by integrating the current density over the strip thickness. The 
conditions under which these assumptions are fulfilled are discussed in Il35l l36l . 
Regarding the magnetic field around a thin superconducting strip, one expects 
that the 2-dependence of the current density can be neglected if the thickness is 
much less than the observation distance, z ^> d. The thin strip approximation is 
therefore expected to be valid for superconducting films with typical thicknesses 
of a few hundred nanometers and trap heights above several microns. We confirm 
this expectation by comparing to numerical calculations that take into account a 
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finite thickness E 

We calculate first the current distribution in the strip and then the outer mag- 
netic field using the Bio-Savart law. In this section, we focus on strips in the 
Meissner state where magnetic flux does not penetrate. Any magnetic field per- 
pendicular to the strip surface is completely shielded by appropriate screening 
currents. If a bias field B is applied parallel to the wide strip surface, there are 
no screening currents in an infinitely thin strip, and the sheet current is given by 
the transport current alone. This situation is changed when the finite thickness is 
taken into account. 

For an infinitely thin strip, the sheet current profile along the a>axis as well 
as the field distribution around the strip can be calculated analytically. The sheet 
current density is given by 05113611 

J(x) = -ML- (1) 



y/w 2 — 



x 



The magnetic field above the strip is presented in Figf2ttop). We actually plot its 
modulus |B(r) | since this provides the magnetic trapping potential in the adiabatic 
approximation. Due to symmetry, the minimum modulus of the total field occurs 
on the 2-axis (above the centre of the strip) where the field caused by the transport 
current / is parallel to the x-axis: 

bm^^ L dx ' = **i (2) 

xK ' ; 2vr 2 J y/ w 2 -x''Hx' 2 + z 2 ) Ww 2 + z 2 

This can be compared to the field above a normally conducting strip. Here, the 
current distribution can be taken as uniform, and the Biot-Savart law gives 

normal wire: BJx = 0, z) = arctan(w / z) (3) 

2tvw 

The field profiles, described by © and © as well as the numerical results for two 
thickness/width ratios are plotted in Fig|3] By adding a homogeneous bias field 
along the x-axis with value B = —B x (x = 0, z t ), a magnetic quadrupole trap 
is formed at height z t . We see that for a superconducting strip, the required bias 
field is smaller than for a normal conductor, by a factor 2/% at small trap height. 
This is a result of the different distribution of current density. Both wires behave 
practically the same above heights z > 2.5 w. Similar results are obtained for a 
cylindrical wire (see Appendix lAj) . 



'These calculations are outlined in Appendix B. For the sake of faster convergence and to avoid 
singular fields, we have 'rounded' the edges of the wires with a radius of curvature r w w/32. 
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Figure 2: (top) Magnetic trap ("side guide") created by a superconducting rectangular 
wire with a transport current, placed in a horizontal bias field (as shown in FigJT]). We plot 
the modulus of the magnetic field in units of hqI / {2ir 2 w). For / = 1 A and w = 1 mm, 
this unit corresponds to (2/tt) G. (top left) infinitely thin wire (analytical calculation); 
(top right) finite thickness d ~ 0.04 x 2w (numerical calculation), (bottom) Cross-sections 
along the z-axis and x-axis through the trap centre (x = is chosen in the middle of the 
top wire face.) Solid lines: analytical results for an infinitely thin strip in the sheet current 
approximation; dashed lines: numerical calculation for a finite thickness. Due to the finite 
thickness, the trap is shifted closer to the wire surface (compared to the case d = the 
trap height is reduced from z t = 0.76 w to 0.64 w). Half of this shift can be explained 
by measuring the distance from the center of the wire, and shifting the dashed curve by 
d/2 would make it coincide with the solid curve at large distances. The bias field is 
Bo = — 2.5 hqI / {2ir 2 w). The numerical calculations use a thickness d « 0.04 x 2w, 
with rounded corners (radius 0.031 w). 
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Figure 3: (left) Magnetic field created by a supercurrent / in a rectangular wire as a 
function of the observation distance z (in units of the strip half-width w, measured from 
the wire top surface). Solid lines give analytical results for d = 0: magenta - normal 
metal strip, green - superconducting strip (superimposed diamonds: numerical calcula- 
tion). Dashed lines (with symbols) give numerical results for superconducting strip of 
different thickness: blue - d = 0.04 x 2w, and red - d = 0.13 x 2w. (right) Bias field 
required to form a trap at distance z% above a strip wire. This plot differs from the left 
one only for wires of finite thickness where screening currents caused by the bias field are 
different on the top and bottom surfaces. For this case, the trap height is slightly reduced. 

Figsl2l[3]also illustrate the impact of the finite thickness of the superconduct- 
ing wire. The bias field induces screening currents on the top and bottom faces 
of the wire and makes the currents on the top and bottom sides differ (Fig{TJ top 
right). The effect is, however, quite small for the parameters that we explored (up 
to thickness d ~ w /A). For example, the magnetic field gradient in the trap centre 
is practically unchanged (Fig©. The main tendency is to bring the trap position 
closer to the wire surface (at fixed bias field), see Fig|2l Conversely, the minimum 
bias field required to create a trap at nonzero height is lowered (Figj3l right). 

3.1.2 Vertical bias field 

The influence of a magnetic field on the current distribution in a flat superconduct- 
ing wire is most pronounced when the field is perpendicular to the wide surface 
of the wire. For this reason, different loading procedures have to be designed for 
superconducting atom chips compared to normal ones. Typically, a chip trap is 
loaded by transporting (positioning) an atom cloud along the x-axis adjusting bias 
fields and currents. In a normal-metal atom chip this is performed using a verti- 
cal bias B z . Above a superconducting wire, the procedure must be altered as the 
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wire builds up significant screening currents to shield its interior from a vertical 
field. The field profile above the chip becomes significantly non-uniform, leading 
to potential barriers that have to be taken into account for the cloud transport. This 
effect is maximal for superconductors in the Meissner state and provides the most 
striking difference to a normal conductor. 

As an example, we illustrate in Figj4]the magnetic field above a superconduct- 
ing strip placed in a vertical bias. Magnetic field maxima occur at the edges of the 
strip (FigHJleft)) that exceed significantly the applied bias field near the wire sur- 
face. This field increase is reduced to about 15% at a height of 1.5 w, but should be 
taken into account for loading the trap. Note that in order to maintain the Meissner 
state, the total magnetic field should be less than the lower critical field B cl in any 
point of the superconductor surface. The vertical magnetic field near the edges 
of the superconducting strip increases proportionally to the ratio width/thickness 
(2w/d) and may become significant for a wide strip ll37l . The same magnetic 
field increase near the strip edges is induced by the transport currently This field 
concentration can result in the partial transition of the sample into the mixed state 
even if the current and magnetic fields are far from the critical values. Magnetic 
traps based on superconducting wires in the mixed states are discussed in the next 
section. 

3.2 Superconducting guiding wire in the mixed state 
3.2.1 Side guide 

Let us consider a side guide trap realized by a type II superconducting wire in the 
form of a strip (Fig[T|), carrying a transport current / in zero external field. To 
calculate the current distribution in the strip we use the Bean critical state model 
[|39l : the current (area) density can only take three different values: ±j c or 0. 
Following Brandt 11351 . the sheet current (integrated over the thickness of the strip) 
is determined as J(x) = (d + — d )j c , where d + and oL are the history-dependent 
thicknesses of the regions carrying +j c or —j c , respectively (d + + gL < d). The 
sheet current cannot be higher than its critical value J c = dj c . This value is 
achieved in regions near the strip edges. In the central part of the strip, a field-free 
region exists that is shielded by the current-carrying domains from the magnetic 
field, similarly to the Meissner state. In this model for a type II superconducting 

2 This is a rough approximation which takes into account only the demagnetisation factor. A 
more accurate analysis is given, for example, in ||38l . 
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Figure 4: Modulus of magnetic field above a thin superconducting strip placed in a 
vertical magnetic field of magnitude Bq. (left) Field at fixed height, |B(x, z = 0.16io)|, 
given in the units of the applied field Bq. (right) Field |B(x = 0, z)\. The position x = 
corresponds to the strip centre above which the field is oriented perpendicular to the strip; 
due to symmetry and continuity, it must go to zero at the wire surface. The lines compare 
calculations for strips of different thickness d. Infinitely thin strip: analytical calculation 
based on Ref.[35]; d > 0: numerical calculations. 



strip, the distribution of the sheet current is given by the expression 113311 : 



(2J c /ir) arctan \ 1% — K, \x\ < b ... 
J(x,J c ,I) = { \ cl 1 V & 2 ~* 2 > ' (4) 

J c , o < \x\ < w 



1 /2 

where b — w (1 — I 2 /I 2 ) , I c = 2wJ c is the maximal value for a supercon- 
ducting current in the strip (critical current); it is achieved at full magnetic field 
penetration, 6 = 0. The distribution of the sheet current is presented in FigJ5](top 
left). Eq.© applies when the transport current / has been increased from zero 
(virgin state) for the first time (cyclic current change is discussed below). 

Due to symmetry, the side guide is located above the strip centre when the 
bias field is parallel to the .x-axis. FigJ5] (bottom left) presents the magnetic field 
produced by a current / on the axis x = of an infinitely thin strip of width 2w 
in the mixed state. The dependence of the magnetic field on the distance z was 
calculated numerically for four ratios of the current to the critical current I / I c 
(the same as in Figj5](top left)). The data should be compared to Figj3] where a 
superconductor in the Meissner state and a normal conductor are considered. The 
magnetic trapping potential that is formed in combination with a horizontal bias 
field (of magnitude 2.5fi I / \2n 2 w)), are shown in Fig J5] (bottom right). These 
plots demonstrate the tendency: the closer the current to the critical current, the 
closer are the trap parameters to the values of a normal metal trap. 
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Figure 5: (top left) Sheet current distribution in a thin strip of type-II superconductor as a 
function of the transport current / which is increased from zero (virgin state). The curves 
are plotted for I / I c = 0.2, 0.5, 0.85, and 0.95 where I c = 2dwj c is the (total) critical 
current; (top right) Magnetic field vs. distance and current above the centre of the strip, 
x = 0, with no external (bias) field. For a critical current density of j c = 10 11 A/m 2 
the field reaches a maximum value of hqI c / i {2'k 2 w) « 127 G{d/1 /im) at I = I c and 
at the wire surface, (bottom left) Magnetic fields calculated as a function of distance 
for different I/I c ratios, as given in the legend: black - Meissner state; red . . . cyan - 
I/I c = 0.2 . . . 0.95; magneta - normal metal. The field is normalized to [IqI /(2tt 2 w), 
as in Fig|2] The calculations were made numerically using the current distributions given 
by ©. (bottom right) The "mixed-state" magnetic potentials as a function of distance 
z at x = are presented for the same I/I c ratios. The horizontal bias field equals 
2.5/i //(2vr 2 ^). 
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It should be noted that a superconductor in the mixed state is essentially a non- 
linear system. As can be seen from the 3D-plot presented in Fig j5] (top right) the 
magnetic field is a nonlinear function of the current, in contrast to both a normal 
conducting wire and a superconducting one in the Meissner state. This happens 
because the shape of the current distribution depends on the ratio I/I c [see Figj5] 
(top left)]. The calculations show that for low currents (/ < 0.2 I c ), the nonlinear- 
ity is negligible and the field distribution around the type II superconducting strip 
may be described by the expressions obtained for a strip in the Meissner state in 
Section 13.11 In the opposite case, / = I c , the current density equals the critical 
value over almost the whole strip width. The magnetic field around the strip can 
then be calculated as for a normal metal, with a spatially uniform current density. 




Figure 6: Illustration of trapped magnetic flux after ramping up and down a transport 
current in a type II superconducting wire in the mixed state. We plot the magnetic field 
produced by a transport current / = 0.85 I c above the centre of a strip wire (curve 1), 
and the magnetic field caused by the flux trapped in the strip, after decreasing the current 
again to zero (curve 2). The magnetic field is given in units of hqI c / (2k 2 w). The inset 
shows the remnant sheet current J(x) in the strip after the current cycle, normalised to 
the critical sheet current J c . The calculation follows Ref. |[35l . 



3.2.2 Trapped magnetic flux and related noise 

The capability of type II superconductors to freeze a magnetic flux results in the 
irreversible behavior of the current profile in the mixed state strip under cyclic 
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changes of the external fields or transport currents. Fig J6] presents the magnetic 
field above the superconducting strip after ramping up the transport current to 
I = 0.85 J c (curve 1) and decreasing it again to zero (curve 2). The inset shows 
the inhomogeneous current distribution that is left in the strip as a result of the 
freezing of the magnetic flux in this current cycle: in the centre, the current keeps 
flowing in the direction of the peak transport current, while near the edges it flows 
backwards. At short distance (z < 0.2 w) the frozen magnetic field is about 15% 
of the field produced by the maximum current. The ratio of these fields decreases 
approximately like 1/z 2 and is about 1% at z = 3w. These calculations were 
made for zero external (bias) field in the framework of the Brandt model P51 . 

The effect of the frozen magnetic field has to be taken into account for the 
proper loading and control of the magnetic trap. The frozen flux can also be the 
main source of magnetic noise generated by the guiding wire because of the mo- 
tion of the trapped vortices. Other effects which can influence the atom cloud 
above mixed state wires are the variations of the critical current density due to lo- 
cal changes in the magnetic field or due to structure inhomogeneities. In addition 
to the vortex lattice, these effects could corrugate the trapping potential at low trap 
heights, z t < w. Changes on slow time scales may occur due to the re-distribution 
of the frozen flux iHOlHTTl. 

The magnetic noise (flux noise) due to the motion of the trapped vortices is 
much higher than the noise in the Meissner state. The flux noise was investigated 
in many experimental p2ll43ll44| and theoretical p2ll44ll45l works. It was shown 
that it is closely related to the mechanism of vortex pinning which in turn depends 
on material properties and preparation technology of the sample. The flux noise 
intensity is determined by external conditions as well. When the temperature, 
current and magnetic field are far from the critical values T c , J c , and B c2 , the 
noise arises from thermally activated, mutually incoherent hoppings among pin- 
ning centres. For this regime of "flux creep", SQUID-measurements near high-T c 
films demonstrate that the noise is strongly dependent on the film quality ll43l . In 
addition, the flux noise picture is more complicated: the noise level depends on 
the magnetic pre-history, it follows not only the number of trapped vortices, but 
also the specific spectrum of metastable trapped states. In particular, a long-lived 
noisy state may occur in a superconducting wire after a pulsed transport current if 
the current is not small compared to the critical value Il43ll48l . 

The influence of the magnetic pre-history is obviously harmful for control- 
ling cold atom traps near a superconducting surface. Further to the additional 
flux noise, the vortices trapped in the superconductor produce magnetic disorder 
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because of the arbitrary locations of the pinning centres. Note that cooling of 
the wire through the superconducting transition in an external magnetic field may 
also lead to the freezing of vortices. SQUID studies of the remnant noise due to 
vortices trapped in a weak field show that the flux noise spectral density at low 
frequencies decreases linearly, as the magnetic field in which the superconducting 
transition occurred, is lowered Il43ll49l . In fact, the remnant noise is proportional 
to the number of vortices trapped in the sample in agreement with the Dutta, Di- 
mon and Horn model ll47l . This number is also determined by the shape of the 
superconductor: the narrower the strip, the less vortices are frozen. As was shown 
in Ref. ll50l . a niobium strip of width 2w = 10 fim placed in a vertical magnetic 
field smaller than ps 0.5 G does not freeze vortices while being cooled through T c . 
At the same time, for a width of 100 /im, the critical field for vortex freezing is 
less than 4mG ll50l . These observations indicate the advantageous nature of nar- 
row superconducting wires for atom chip experiments which require particularly 
long lifetimes and coherence times. 

To summarise the possibility of using mixed state wires for atom guiding and 
trapping, the advantage over the Meissner state is the higher current and thus 
the tighter confinement, while the disadvantage is the expected static potential 
corrugation and higher noise. It may even turn out that their performances are 
comparable to other, normally conducting materials that have been suggested re- 
cently llFTl 1521 . but more detailed investigations are needed here. On the other 
hand, it may be possible to trap cold atoms close enough and long enough so that 
they may be considered as a probe of the mixed state of type II superconductors. 
The low-frequency flux noise in the mixed state could be detected by the spin 
dephasing rate of cold atoms, as was noted in ll53l . The high sensitivity of cold 
atoms to a disturbance of the magnetic potential could be used for visualization of 
the static disorder produced by frozen vortices, analogously to the current static 
scattering in normal atom chip wires ll54l[T8l[T9ll . Atom clouds having high spatial 
resolution (3 /im) combined with excellent sensitivity to magnetic field (4 nT) 11541 
(or even better than that Il55l0 could provide complementary information about the 
distribution and dynamics of vortices. According to the review in Ref. ll56ll . such a 
combination of resolution and sensitivity would be one of the best among various 
vortex observation methods. 
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4 Comparison between normal and superconduct- 
ing magnetic traps 

Our calculations show that magnetic traps can be created in atom chips with su- 
perconducting wires. The main differences to normal wires are the inhomogeneity 
of the current distribution and nonzero screening currents induced by bias fields. 
The current density in a superconducting wire is smaller in its centre compared 
to a normal conductor (at the same total current). For this reason, the magnetic 
field near the wire surface is weaker, and a side guide trap (in a given parallel bias 
field) is closer to the surface (see Figj3]). 

Let us now analyse in more detail the difference in trap parameters between 
side-guide traps created by the normal and superconducting guiding wires. We 
also compare the Meissner and mixed states with respect to their trapping "capa- 
bilities". Note that we do not take into account here the bending of the wires into 
"U"- and "Z"-shapes for 3D traps (see, for example, HI El El). Calculation of the 
current density in bent superconducting wires is complicated because it requires 
to solve a three-dimensional problem even in a planar configuration [|20l . Our re- 
sults can be applied to the central part of the wires, sufficiently far away from the 
bends. 

Two figures of merit describing the confinement of cold atoms in the magnetic 
trap are used for this comparison: the magnetic field gradient at the trap centre and 
the depth of the trapping potential. The trap depth is determined as the minimal 
height of the (total) potential barrier, from the trap centre to either the surface or 
away from it. Here, the gravitational potential is taken into account. We adjust 
the trap height by setting the bias field to the required value and take the same 
geometry and total current for a fair comparison between superconducting and 
metallic wires. 

Results for typical flat wires are shown in FigJTl Comparing superconducting 
and normal guiding wires, we see that the field gradients (left scale) and the trap 
depths are essentially the same for trap distances z t > 2 w, similar to what is seen 
in Figf3] The difference is maximal in the range of low trap heights, where one 
becomes sensitive to the weaker magnetic field and current density in the super- 
conductor. This result is in qualitative agreement with the numerical simulations 
of confinement parameters of specific superconducting traps performed in ll2~0l . 
The finite thickness of the wire does not significantly change the picture, as can 
be seen from the numerical data plotted in FigjTfb) (symbols). For example, an 
infinitely thin superconducting strip creates a magnetic field gradient along the 
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Figure 7: Comparison of trap parameters (field gradient and trap depth) calculated for 
side-guide traps with superconducting and normal metal wires, (a) Superconductor in 
the mixed state (YE^C^O/-^ at 4.2 K) and normal metal (black and red lines, respec- 
tively). Calculations were performed analytically for infinitely thin guiding wires of width 
2w = 10 fim, and current I = 200 mA w0.1J c (assuming a thickness d = 300 nm). The 
trap height is varied by adjusting the bias field. Dashed lines: absolute value of mag- 
netic field gradient \dB/dz\; solid lines: trap depth, (b) Superconductor in the Meissner 
state (niobium at 4.2 K). Lines: analytical calculation for an infinitely thin wire of width 
2w = 10 /an, carrying 20 mA of current. Black and red squares: numerically determined 
trap parameters for a strip with thickness d = 2 fim. The numerical calculation takes 
into account the additional screening currents created by the bias field (negligible for an 
infinitely thin wire). Dash-dotted (horizontal) lines: minimum field gradient required to 
stabilize the trap against gravity and minimum trap depth to confine a cold atomic sample 
at 1 fiK. The trap depth is calculated from the height of the potential barriers towards or 
away from the chip, including gravity ("upside down" setup with the trap below the chip). 
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z-direction at the trap centre (height z t = 0.75 w) that is about 1.8 times less than 
in a normal metal trap. Along the x-axis the ratio of the gradients (normal metal 
to superconductor) is about 1.3. 

The data presented in Fig|7] are calculated for the high-T c superconductor 
YBa 2 Cu 3 7 „ ( 5 (YBCO) and for niobium (Nb) in the Meissner state. If the YBCO 
film is used at 4.2 K, the assumed current of 200 mA is about 10% of the critical 
value (critical current density j c ps 7 ■ 10 11 A/m 2 ). Nevertheless, we expect the 
wire to be in the mixed state, at least at its edges where the current density exceeds 
the values permitting the Meissner state. For niobium, the critical current density 
is lower, and this is why we reduce the transport current to 20 mA. The horizontal 
lines in Fig j7] mark typical criteria for reliable trapping of atoms at a temperature 
of 1 yuK: the trap depth should exceed 10 /jK and the gradient should be high 
enough to protect the atoms from gravity's pull (corresponding to 15.3 G/cm for 
87 Rb atoms in the \F = 2, nip = 2) state, where F is the total spin and mp its 
projection on the local magnetic field). 

According to Fig|7£a), we predict that an atom chip based on a YBCO super- 
conducting strip can trap cold atoms in a wide range of trap heights 0.2 4- 300 fim. 
For niobium, the range is smaller (1 -h 75 fxm) due to the lower guiding current 
(Figl2b)). The trap parameters are still high enough, however, to successfully trap 
cold atoms, both in the Bose-Einstein condensed phase and above. It is also seen 
in Fig j7] that the trap parameters in the closest vicinity of the surface (z < 5 /im) 
are worse than for a normal strip. However, a significant gain in the lifetime due 
to the reduction of magnetic noise near the superconductor makes this trap design 
more attractive. 

5 Conclusion 

We have presented a theoretical analysis of superconducting atom chips. Our 
methods have been both analytical and numerical; they complement the results 
recently reported in Ref.[20]. In particular, the analytical expressions given here 
can be used for semi-quantitative estimates and to identify scaling laws for chip 
design. We confirm the possibility of trapping cold atoms in a wide range of 
distances (0.2 300 fim) with the same wire. This analysis takes into account the 
specific behavior of superconductors that carry a transport current in a magnetic 
field. These peculiarities enforce modifications in the loading procedure and the 
control of the atom cloud. The application of superconductors to atom chips may 
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enable improved atom optics with suppressed effects of noise, as well as novel 
insight regarding the noise and current distribution in superconductors. 
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A Side guide magnetic trap based on cylindrical su- 
perconducting wire 

Let us consider a conductor in the form of a cylinder (radius R) with a DC trans- 
port current / in an external magnetic field B (FigJH]). A normal metal wire does 
not influence the external DC field so that the total field is the superposition of the 
homogeneous bias field and the field produced by the current. The total field is 
zero at 9 = and at the trap height 

z t = r t -R = -^j- - R (5) 
2ttB 

A magnetic field cannot penetrate a superconductor in the Meissner state. If 
the radius is much larger than the London penetration depth, the magnetic field 
is zero inside the cylinder. The magnetic field around the superconductor without 
transport current can be described by the Laplace equation for the scalar magnetic 
potential ip. The boundary conditions are: the radial component of the magnetic 
field equals zero at r = R; at r — > oo, the magnetic field coincides with the 
external one. The solution is 

B r = -d r ip = B ( 1 - ) sin 9 (6) 




(7) 



where B r and B e are the radial and azimuthal components of the magnetic field. 
The magnetic field configuration around the superconductor thus differs from a 
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Figure 8: Cylindrical wire and coordinate system. 

normal metal even if a transport current is absent. Due to the Meissner effect, the 
magnetic field near the superconductor at 9 = is increased. 

If the cylinder carries a transport current, the outer field is the same as for a 
normal conductor, due to the cylindrical symmetry. We find that the minimum of 
the total field (Fig|9£left)) is located at 9 = and distance 

„ 1 + y/1 ~4:h 2 „ , 2ttRB 

z t = r t -R = R W — R, h= -A (8) 

in fi l 

from the wire surface. Because of the screening of the bias field by a supercon- 
ductor, the same trap height in a superconducting chip is achieved with a lower 
bias field than in a normal conducting chip (see FigHright) and also FigOJright) 
for a rectangular wire). This reduction reaches 50% at small trap heights. The 
difference decreases with an increase of this height and practically vanishes at 
zt > SR. 

B Numerical calculation of magnetic field around a 
finite thickness superconductor 

We consider a wire that is infinitely long in the y-direction and ignore here bound- 
ary effects (these become relevant for "U"- and "Z"-shaped wires, of course). By 
symmetry, the magnetic field is independent of y and lies in the icz-plane (see 
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Figure 9: (left) Modulus of the total magnetic field around a superconducting cylinder 
with a transport current / in a bias magnetic field Bq = 0.4 x hqI /2ttR. The field 
minimum is located at the point 9 = 0, z t = 0.8 R (r = 1.8 R). (right) Required bias 
field as a function of the trap height for a cylindrical wire, based on Eqs.© and ([8]): solid 
line - superconductor; dashed line - normal metal. In both graphs, the magnetic field is 
normalised to /j, I/2ttR = 2G(//lA)/(E/lmm). 

Figffl). Inside the superconductor, the field is zero. Outside, we describe the field 
in terms of a scalar potential ip and a vector potential A = e y A 

B(x) = - W(x) + V x A(x) . (9) 

We shall see below that the scalar potential describes the field caused by an ex- 
ternal magnetic field (as in Appendix A), while the vector potential gives the field 
caused by a transport current in the superconducting wire. The introduction of 
two potentials may seem superfluous, because one could work, outside the wire, 
with the scalar potential only. But tp would then become a multivalued function 
for a nonzero transport current. 

The magnetic field normal to the superconductor surface is continuous and 
must therefore vanish. (This corresponds to the radial derivative of ip in Appendix 
A.) Writing n(x s ) for the normal at surface point x s , we have 

n(x s ) ■ B(x s ) = -n(x s ) ■ V^(x s ) - (n(x s ) x e y ) ■ VA(x s ) 
^(x s ) dAjx,) Q 
<9n(x s ) dt(x s ) 

where d/dn(x s ) and d/dt(x. s ) are the normal and tangential derivatives at x s . 
(The local tangent vector is t = nxe r ). We construct the potentials if) and A such 
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that both terms d^i s jdn and dA s /dt in Eq.dTOb are zero. Note that dA s /dt = 
indicates that A(x s ) is a constant on the surface which we denote by A . The 
tangential magnetic field is nonzero at the (outer) surface of the superconductor; 
the interior of the wire is screened from this field by a surface current (along the 
?/-axis) of magnitude t • B//i - 

Considering that outside the wire, both divergence and curl of the magnetic 
field are zero, we get 

V 2 ^(r) = 0, V 2 A(r)=0 (11) 

where V 2 is the two-dimensional Laplace operator. These Laplace equations can 
be solved with the help of the Green function 

G(x,r) = - 7 llogr, (12) 

where r = |r — x|, given the values of the potentials and their derivatives on the 
wire surface, and an asymptotic condition at large distance. 

The two potentials behave differently with respect to the current / transported 
by the superconductor. It can be seen by recalling the Ampere-Maxwell law (in 
the static limit) 

^ / = /B(x)-da=-jf ^da, (13) 

where / is the total current flowing through the superconductor, S a closed curve 
including the cross-section of the superconductor (with oriented line element da), 
and d/dn is the derivative normal to the curve S and pointing 'outside'. Hence, 
the vector potential A is proportional to the transport current. Conversely, the 
scalar potential asymptotically goes over into ip — > Vext( r ) = — r ■ B where B is 
the homogeneous bias field. This asymptotic condition forces ij) to be proportional 
to B . 

We solve the Laplace equations (PTTI) for the potentials M = ip, A in terms of 
a surface integral equation (r outside the wire) Il57ll58l 

M(r) = «_(,, - / s d a( x) (G(x, r) ^W - ^M(x)) , (14) 

where M cxt is the external potential (nonzero only for ip), S is now the circum- 
ference of the wire (with scalar line element da(x)), and the normal derivative 
d/dn(x) points 'outside' the wire. For M = ip, we have dip(x)/dn(x) = on the 
surface according to Eq.dTOb; then there is only one unknown, ^>(x), in Eq. (fT4)) . 
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When M = A, we can set under the integral A(x) = A , constant on the surface, 
also according to Eq.dTOb. Then the integral of the second term in the bracket can 
be shown to vanish, and there is only one unknown left, dA{x)/dn{x), which is 
actually the surface current density, see Eq.(fT3l. We always choose A = 1 in the 
beginning, compute the total current by Eq.(TT3l (giving a coefficient like the wire 
inductance) and re-scale A(x) and dA/dn(x) to get the desired current. 

The potential M and its derivative on the surface are obtained by letting r — > 
x s in Eq. (fl4|) approach the surface. Note that we touch the singularity of the 
Green function G(x, x s ) and its normal derivative dG(x, x s )/dn(x) under the 
integral when x = x s . To handle this, we expand the integrand within a small 
neighborhood s of length da around x s and perform the integration, giving 

/da(x) ^M( X ,,(i + ^)M( X< ). (1 6) 

Here the angle d(p(x s ) involves the radius of curvature i?(x s ) of the surface: 
d0(x s ) = da/R(x s ), it describes the angle subtended by the surface element 
as seen from the centre of curvature. The logarithmic correction of Eq.(fl~5l) and 
the curvature correction (PT6l) greatly improve the convergence of the numerical 
calculations. 

The other parts of the surface integral, excluding the point x s , are discretized 
in the usual way, mapping the integral equation into a linear system ||59ll60ll . For 
example, a rectangle 20.64 /zmx 0.84 /im with rounded corners (R = 0.32/im) 
is typically discretized into 420 surface elements along its circumference. This 
leads to a linear system with 420 2 matrix elements. Once we have the potential or 
its derivative on the wire surface, the field outside the wire is found from Eq. (fT4)) . 
We have checked the convergence of the numerics, for different discretizations of 
the wire surface, against the exact solution of Appendix A for cylindrical wires 
and against the solution of Ref . 11611 for rectangular wires with sharp corners. 
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